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Abstract 



We propose a new macroscopic model derived from the classical nonlinear Boltz- 
mann equation. A set of partial differential equations is obtained easily. The 5N 
unknowns Naj{t,x) depend on the time and space coordinates, a = 1,2,.... A and 
j = 0, 1.., 4. The index a labels a subset Cq of a partition of the velocity space and 
the five indexes j correspond to the collision invariants. The unknowns are related 
to the distribution function /(t,x, solution of the Boltzmann equation, since 



This new model guarantees the conservation of the mass, momentum and energy. 
We prove that the set of equations coincides with the set obtained applying the 
discontinuous Galerkin method to the Boltzmamn equation [1]. 




Aa,i(t,x) 
Aa,2(t,x) 

A^a,3(i,x) 




MSC-class: 76P, 82C40 (Primary) 65M60 (Secondary) 



1 



1 Introduction and basic equations 



We consider the classical nonlinear Boltzmann equation [2], [3] for neutral monatomic 
gases 

| + «-|-Q(/./)^ (1) 

The one-particle distribution function / depends on time t, position x and velocity ^. 
We denote by X C M'^ the domain of the spatial coordinates x. The velocity space is M^. 
The collision operator is defined as follows 

Q{fJ)= [ I I w{tiAi'.OU'fl-fh)di.di'dC (2) 

JR3 J]r3 J^i 

Here, as in the following, to simplify the notation, often we omit to write the variables t 
and X, explicitly. Moreover, it is here understood that 

/' = /(t,x,o, /: = /(t,x,o, / = /(t,x,^), /, = /(t,x,^,). 

The kernel W of the collision operator is defined by 

w{i, o = Kin . V, I vi) + - r - a m' + - I' - (3) 

where 

" = 1^11 V = ^-e. (4) 

The function K is related to the interaction law between colliding particles. The Dirac 
distributions guarantee momentum and energy conservation during the binary collisions. 
The collision operator Q{f, f) of the Boltzmann equation ([T]) is usually written in different 
way, since Dirac distributions are used to reduce the collision integral to a five-fold integral. 

Let : — i> M be a measurable function. Now, we assume that generic test functions 
(p depend only on the velocity If we multiply both sides of the Boltzmann equation ([1]) 
by a test function 0(^) and we integrate with respect to the velocity ^, then we obtain 
the equation 



d_ 

dt 



[ /(t,x,O0(Orf^ + Vx/ /(t,x,^)C0(O^^= / Q{fJ){t,^:Omd^. (5) 

jRf Jr3 Jr3 

For any function (f), we obtain an equation of type ([5]). This procedure recalls the classical 
moment method, but, in our framework, the test functions will have compact support. 

We recall a well-known result. Assuming the existence of the integrals, the right hand 
side of Eq. ([5]) can be written as follows 

If f f [ W{tU^',0[<t>' + €-<l^-<P*]ff*d^'dej^,d^. (6) 

^ JR3 Jr3 J^3 
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Now, if we define 

W;^,e)=/ / w{tu^',am')+HO-m-m)]d^'dc (r) 

then the integral becomes 

/ Q{fJ)md^ = l f [ /C(0;^,e)/(O/(^*)^l*rf^ (8) 
and Eq. writes 

^/ f{t,^,^)md^+vj f{t,^,^)^md^ 

= \ I f n4>;ti*)f{t.^,i)fit.^.i*)di.di. (9) 

We note tliat the function lC{(t)\ ^, ^^,), which is related to the function 0, depends only on 
the (^,^*) variables and plays the role of a kernel of the integral operator (|H]). We recall 
a simple but important result. If ip denotes one of the collision invariants 1, ^, then 
we have 

/C(V';^,e) = o VC^gM^. (10) 

We denote by {tfjj : j = 0, 1.., 4} the ordered set {1, ^1,^2, ^3, I^P}; where ^1, ^2, ^3 are the 
three components of the vector ^. 



2 The macroscopic model 

As in the numerical treatment of a kinetic equation by means of finite differences or 
elements, we requires a bounded domain for the velocity space. To this scope, we introduce 
a suitable characteristic function in the kernel of the collision operator, such that there 
exists a bounded domain f2 C so that, if /(0,x, ^) = for every ^ ^ and x G X, 
then fit, X, ^) = for every ^ ^ f2, x G X and for all time t. Let be a suitable positive 
real number. We define the function : x — M as follows 



and 



1 if < 
otherwise 



WE{t eir , o = xe{^, e) Wit t\^', o ■ (11) 

The modified kernel We has the same fundamental properties of the true kernel W and 
guarantees that if two particles have velocities ^ and t before the impact, such that 
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+ < E, then the velocities ^' and ^' of the particles, after the impact, will satisfy 
the inequality + < E. We can choose 1] = G : < E}, and we consider 
the Boltzmann equation with the modified kernel We- 

Now, we choose measurable subsets Ca {a = 1,2,,.., A^) of Q such that 



N 



C^cn Va , n = Va ^ /? and [j = ^ ■ 



a=l 



We remark that we have not introduced any constrain on the size and shape of the cells; 
so, there is a great arbitrariness in the decomposition of the domain Q. We denote by Xa 
the characteristic function on the set C^. 

We consider the 5N test functions 



^(l)Xa(O for iji^) = 1,6,6,6, III' and a = 1,2, ....N . 
Therefore, Eq. furnishes the 5A^ equations 



(12) 



(13) 



Q, Jn 



where a = 1, 2, ....A^ and j = 0, 1.., 4. 
We choose the 5A^ scalar integrals 

A^„,,(t,x) 



/(t,x,OV^,(0 d^ 



(14) 



as the only independent variables and we look for a reasonable closure of the system (IT^ . 
The physical meaning of the new variables is evident; for instance, Na.o{t, x) is the density 
of particles at time t and position x having velocity belonging to the set Ca- 

A very simple recipe for a closure is the following. 
Firstly, we treat the drift term of Eq. f lT^ . We consider a generic cell Ca and we look for 
an approximation in Ca of the function ^ ipj{$,) {j = 0, 1.., 4) of this type 



where a, 



(a) 
kj 



k=0 

{a) 
^kj 



i^ ^.(1) ~ , Mi) for every ^ G and ^ = 1, 2, 3 , (15) 

is a constant array to be determined. 



J 2 



(q) 
^kj 



4 



Now, we consider the collision term of Eq. f lT^ . Since 
^ Jn Jn 



N N 



8=1 7=1 •^'^13 •^'^1 



13=1 7 

we look for an approximation of this type 

4 4 



XE{i. e) X.\ t ^ E E tknM^) ^ni^*) e^e^y e X . (16) 



k=0 n=0 



Here, fe^-^f^'' are the numerical parameters to be determined. 

If we assume reasonable the approximations given by Eqs. ( I15p -f ll6p . we obtain the fol- 
lowing closed system of partial differential equations 

a AT ^ N N 4 i 

^ + E ■ v.iv«. = i E E E E ^..» ■ (17) 

A:=0 /3=1 7=1 k=0 n=0 

At this point, we must give a meaning to the approximations given by Eqs. f|T5|) - f|T6|) . 
For each cell Ca we introduce a vector function r]a{$,)- The components of the five 
dimensional array T7a(^) are functions, denoted by ?7a,j(^); which are linear combination 
of the collision invariants and such that 

/ VaA^)i^M)d^ = S,, (18) 
for every i and j and for each cell Ca- 

We determine the numerical parameters a^"^ and h^^^^\ for every (a, /3, 7 = 1, N) and 
(j, A;, n = 0, 4), assuming that the following equations hold 

/ n.M)iUi)di = Y.^tn VaA^)M^)d^, (19) 

?7/3,i(0 ?77,m(^*) ^(^i Xa, ^, ^*) <i^* C?^ 



4 4 

= E E ^%n^ / / M^) ^n(e) • (20) 

fc=0 n=0 JCp JC~, 

Taking into account Eq. (fT8|) . we obtain 

aSf = / r/,.(0^^.(0^^^, (21) 

= / / ^A^(^) ^7,-(e) Xi.(C, C*) /C(V^, t e) c?^. ciC • (22) 



We note that the parameters a^^ and are numerical constants, which depend only 

on the domain decomposition and the scattering kernel K. Therefore, they do not depend 
on the solutions of the Bolzmann equation. 

It is immediate to verify that, with this choice of the parameters, the macroscopic equa- 
tions ( TTTIl coincide with the equations (Ref. pQ) obtained applying the discontinuous 
Galerkin method to the Boltzmann equation ([1]). Moreover, as a consequence of a re- 
sult of paper [1], it is guaranteed the conservation of mass, momentum and energy for 
homogeneous solutions. 
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